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Abstract 



We report on the first measurement of inclusive ]/\ff elliptic flow, V2, in heavy-ion collisions at the 
LHC. The measurement is performed with the ALICE detector in Pb-Pb collisions at ■JSfw = 2.76 
TeV in the rapidity range 2.5 < y < 4.0. The dependence of the ]/\f/ V2 on the collision centrality 
and on the ]/y/ transverse momentum is studied in the range < pr < 10 GeV/c. For semi-central 
Pb-Pb collisions at ^/snn = 2-76 TeV, an indication of non-zero V2 is observed with a maximum 
value of v'2 = 0.1 16 ± 0.046(stat.) ± 0.029(syst.) for in the transverse momentum range 2 < 
, pj < 4 GeV/c. The elliptic flow measurement complements the previously reported ALICE results 

on the inclusive J/y/ nuclear modification factor and favors the scenario of a significant fraction of 
^ , J/y/ production from charm quarks in a deconfined partonic phase. 



*See Appendix lAl for the list of collaboration members 
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The aim of ultra-relativistic heavy nuclei collisions is the study of nuclear matter at high temperature and 
pressure where Quantum Chromodynamics predicts the existence of a deconfined state of partonic mat- 
ter, the Quark-Gluon Plasma (QGP). Heavy quarks are expected to be produced in the primary partonic 
scatterings and to interact with this partonic medium making them ideal probes of the QGP The measure- 
ment of quarkonium states and hadrons with open heavy flavor is therefore expected to provide essential 
information on the properties of the strongly-interacting system formed in the early stages of heavy-ion 
collisions |Q]]. According to the color-screening model quarkonium states will be suppressed in the 
medium with different dissociation probabilities for the various states. Recently, the CMS Collaboration 
at the Large Hadron Collider (LHC) reported about the observation of the sequential suppression in the T 
sector [3]. The ALICE Collaboration published the inclusive J/y nuclear modification factor Raa down 
to zero transverse momentum (pj) at forward rapidity in Pb-Pb collisions at = 2.76 TeV [4]. The 

Raa compares the yields in Pb-Pb to those in pp collisions scaled by the number of binary nucleon- 
nucleon collisions. The inclusive J/y nuclear modification factor reported is larger than that measured at 
the SPS [5] and at RHIC [ 6, 7Q for central collisions and does not exhibit a significant centrality depen- 
dence. Complementarily, the CMS Collaboration measures the high pi (6.5 < pi < 30 GeV/c) prompt 
J/y Raa in Pb-Pb collisions at ^/s^ = 2.76 TeV in the rapidity range \y\ < 2.4 [8]. The data of CMS 
shows that high pj J/y are found to be more suppressed than low pi J/y and that this suppression does 
exhibit a strong centrality dependence. 

The weak centrality dependence of the J/y Raa at low transverse momentum can be qualitatively under- 
stood with models including full 10, [13] or partial [H, 12] regeneration of J/i/Afrom deconfined charm 



quarks in the medium. The J/y regeneration mechanism was first proposed by the Statistical Hadroniza- 
tion Model (SHM), which assumes deconfinement and thermal equilibrium of the bulk of cc pairs to 
produce J/y at the phase boundary by statistical hadronization only |9|] . Later, the transport models 
proposed a dynamical competition between the J/y suppression by the QGP and the regeneration mech- 
anism, which enables them to describe also the pj dependence of the J/y Raa 12]. These models 
have in common the assumption of deconfinement and some degree of charm quark thermalization. More 
differential studies, like the J/y elliptic flow, could help to assess the charm quark thermalization in the 
medium. 

The azimuthal distribution of particles in the plane perpendicular to the beam direction is an experimental 
observable also sensitive to the dynamics of the early stages of heavy-ion collisions. When nuclei collide 
at finite impact parameter (non-central collisions), the geometrical overlap region and, therefore, the 
initial matter distribution is anisotropic (almond-shaped). If the matter is strongly interacting, this spatial 
asymmetry is converted via multiple collisions into an anisotropic momentum distribution [| 1 311 . The 
second coefficient of the Fourier expansion describing the final state particle azimuthal distribution with 
respect to the reaction plane, v%, is called elliptic flow. The reaction plane is defined by the beam axis 
and the impact parameter vector of the colliding nuclei. 



Within the transport model scenario Hill I12H the observed J/y have two origins. First, primordial J/y , 
which are produced in the initial hard scatterings, traverse and interact with the created medium. During 
this process they may be dissociated. Second, J/y could be regenerated from deconfined charm quarks 
and anti-quarks in the QGP. Primordial J/y emitted in-plane traverse a shorter path through the medium 
than those emitted out-of-plane resulting in a small azimuthal anisotropy for the surviving J/y. Regen- 
erated J/y inherit the elliptic flow of the charm quarks in the QGP. If charm quarks do thermalize in the 
QGP then J/y formed there can exhibit a significant elliptic flow. 

At RHIC energies, the (preliminary) measuremen ts by the (PHENIX) STAR Collaboration of the J/y elliptic 



flow in Au-Au collisions at ,Ann = 200 GeV Itl4l ll5ll are consistent with zero albeit with large uncertain- 



ties in the p T and centrality ranges (0-5 GeV/c) 2-10 GeV/c and (20%-60%) 10%-50%. The measure- 
ment of quarkonium elliptic flow is especially promising at the LHC where the high energy density of 
the medium and the large number of cc pairs produced in Pb-Pb collisions should favor the development 
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of flow and the regeneration mechanisms. 

In this Letter, we report ALICE results on inclusive J/y/ elliptic flow in Pb-Pb collisions at 
TeV at forward rapidity, measured via the }l + pL~ decay channel. The J/y elliptic flow is presented as a 
function of transverse momentum and collision centrality. 

The ALICE detector is described in [16]. At forward rapidity (2.5 < y < 4) theproduction of quarko- 
nium states is measured via the decay channel jU + jU~ in the muon spectrometer [j down to pj = 0. The 
spectrometer consists of a ten interaction length thick absorber stopping the hadrons in front of five track- 
ing stations comprising two planes of cathode pad chambers each, with the third station inside a dipole 
magnet delivering a 3 Tm field integral. The tracking apparatus is completed by a triggering system 
made of four planes of resistive plate chambers downstream of a 1.2 m thick iron wall, which absorbs 
secondary hadrons escaping from the front absorber and low momentum muons coming mainly from % 
and K decays. In addition, the silicon pixel detector (SPD) and scintillator arrays (VZERO) were used 
in this analysis. The VZERO counters consist of two arrays of 32 scintillator sectors each distributed in 
four rings covering 2.8 < r\ < 5.1 (VZERO-A) and -3.7 < 77 < —1.7 (VZERO-C). The SPD, used to 
determine the location of the interaction point, consists of two cylindrical layers covering \r\ \ < 2.0 and 
\r\ \ < 1 .4 for the inner and outer layers, respectively. All of these detectors have full azimuthal coverage. 
The data sample used for this analysis, collected in 2011, amounts to 17 x 10 6 dimuon unlike sign (MU) 
triggered Pb-Pb collisions and corresponds to an integrated luminosity J2fj nt m 70 /lb -1 . In addition to a 
minimum bias (MB) trigger, the MU trigger requires at least a pair of opposite sign track segments, each 
with a pt above the threshold of the on-line trigger algorithm. The pj threshold of the trigger algorithm 
was set to provide 50% efficiency for muon tracks with pj = 1 GeV/c. The MB trigger requires a signal 
in VZERO-A and a signal in VZERO-C. The beam induced background was further reduced offline using 
the VZERO and the zero degree calorimeter (ZDC) timing information. The contribution from electro- 



magnetic processes was removed by requiring a minimum energy deposited in the neutron ZDCs 1 17]. 
The centrality determination is based on a fit of the VZERO amplitude distribution as described in 111 811 . 
A cut corresponding to the most central 90% of the nuclear cross section was applied; for these events 
the MB trigger is fully efficient and the contribution from electromagnetic processes is negligible. 

J/Y candidates are formed by combining pairs of opposite-sign (OS) tracks reconstructed in the geomet- 
rical acceptance of the muon spectrometer. To improve the muon identification, the reconstructed tracks 
in the muon tracking chambers are required to match a track segment in the muon trigger system above 
the pj threshold of the on-line trigger algorithm. 

The J/y/" V2 is calculated using event plane (EP) based methods. The azimuthal angle *P of the second 



harmonic event plane is used as an estimate of the reaction plane angle 111 911 . *P is determined from 
the azimuthal distribution of the VZERO amplitude. The VZERO-C has a common acceptance region 
with the muon spectrometer. Therefore, only the VZERO-A was used for the event plane determination 
to avoid autocorrelations. The J/1//V2 results were obtained determining V2 = (cos2(0 — *P)) versus 
invariant mass (ra w ) lEoTl . where is the dimuon azimuthal angle. In this method, V2 of the OS dimuons 
is calculated as a function of and then the resulting V2(m^ l ) distribution is fitted using: 

V2K/i) = V2 g cc{m m ) + vf g (m^)[l - a(m MjU )], (1) 

where v^ lg and Vj kg correspond to the V2 of the J/y signal and of the background, respectively. v^ kg was 
parametrized using a second order polynomial (see Fig. [T|(b)). Here, Of(m w ) = S/(S + B) is the ratio 
of the signal over the sum of the signal plus background of the distributions. It is extracted from 
fits to the OS invariant mass distribution (see Fig. [T|(a)) in each pj and centrality class. The OS dimuon 
invariant mass distribution was fitted with a Crystal Ball (CB) function to reproduce the J/y line shape, 



'in the ALICE reference frame, the muon spectrometer covers a negative 77 range and consequently a negative y range. We 
have chosen to present our results with a positive y notation. 
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Fig. 1: (Color online) Invariant mass spectrum (a) and (cos2(0 — ¥)) (b) as a function of the invariant mass of 
pairs (fitted with Eq. Q) with 2 < p T < 4 GeV/c and 2.5 < y < 4 in the 20%-40% semi-central Pb-Pb 
collisions. 

and either a third order polynomial or a gaussian with a width linearly varying with mass to describe 
the underlying continuum. The CB function connects a Gaussian core with a power-law tail [21] at low 
mass to account for energy loss fluctuations and radiative decays. An extended CB function with an 
additional power-law tail at high mass, to account for alignment and calibration biases, was also used. 
The combination of several CB and underlying continuum parametrizations described before were tested 
to assess the signal and the related systematic uncertainties. The J/i//" V2 and its statistical uncertainty 
in each pj and centrality class were determined as the average of the Vj 8 obtained by fitting V2(m/m) 
using Eq. [Qwith the various a(m^), while the corresponding systematic uncertainties were defined as 
the RMS of these results. Figured] shows typical fits of the OS invariant mass distribution (a) and of the 
(cos2(c/> — *F)) as a function of ra^ (b) in the 20%^10% centrality class. The systematic uncertainty 
related to the unknown shape of the v 2 (m^) was evaluated by repeating the procedure above using 
either a first order polynomial or its inverse as v^ g parametrization. The largest deviation of the results 
obtained with the three different v^ g parametrizations was conservatively adopted as the systematic 
uncertainty. A similar method is used to extract the uncorrected (for detector acceptance and efficiency) 
average transverse momentum ((pt)™ 001 ) of the reconstructed J/i/a in each centrality and pj class. The 
(px) uncor is used to locate the data points when plotted as a function of transverse momentum. Consistent 
V2 values were obtained using an alternative method 01911 in which the J/i/^raw yield is extracted, as 
described before, in bins of (c/> — *F) and the V2 values are evaluated using a fit to the data with the 
function g^rsm = A [ 1 + 2v2 cos 2 ( <j> — *P) ] , where A is a normalization constant. 

The finite resolution in the event plane determination smears out the azimuthal distributions and leads to 
a lower value for the measured anisotropy 111911 . The VZERO-A event plane resolution as a function of 
the centrality was determined using MB events and the 3 sub-event method H19TI . To estimate the system- 
atic uncertainty from the event plane determination two sets of 3 sub-events were used: first, VZERO-A, 
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Table 1: VZERO-A event plane resolution for the centrality classes expressed in percentages of the nuclear cross 

section 01. 

Centrality (N par t} EP resolution ± (stat.) ± (syst.) 



5%-20% 
20%^0% 
40%-60% 
60%-90% 
20%-60% 



^CNJ 0.3 



283 ± 4 0.548 ± 0.003 ± 0.009 

157 ± 3 0.610 ± 0.002 ± 0.008 

69 ±2 0.451 ±0.003 ±0.008 

15 ±1 0.185 ±0.005 ±0.013 

113 ± 3 0.576 ± 0.002 ± 0.008 
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Fig. 2: (Color online) Inclusive J/y v%{j)j) f° r semi-central (2Q%-40%) Pb-Pb collisions at -^/snn = 2.76 TeV. 
The V2 was measured in the pj ranges: 0-2, 2-4, 4-6 and 6-10 GeV/c and the points are located at the measured 
(px) uncor . 



VZERO-C and the Time Projection Chamber (TPC), with pseudo-rapidity gaps At]voa-tpc=1.9 and 
A?7tpc-voc=0.8; second, VZERO-A, VZERO-C-lst ring and VZERO-C-4th ring, with pseudo-rapidity 
gaps Atjvoa-voci=4.5 and At7voci-voC4=1-1- The differences between the event plane resolution for 
VZERO-A obtained from these two sets of sub-events are taken as systematic uncertainties. Since V2 
is measured here in a wide centrality class, the resolution must reflect the distribution of events with 
a J/y within the class. Therefore, the event plane resolution for each wide class was calculated as the 
average of the values obtained in finer centrality classes weighted by the number of reconstructed J/y. 
Table [TJ shows the corresponding resolution for each centrality class which is applied to the results re- 
ported in this Letter. 

The J/y reconstruction efficiency depends on the detector occupancy which could result in a bias of the 
V2 measurement. This effect was evaluated by embedding azimuthally isotropic simulated J/y — >• }X + pL~ 
decays into real events. The measured V2 of those embedded J/y was found not to deviate from zero 
by more than 0.015 in all the centrality and pj classes considered in this Letter. This value is used as a 
conservative systematic uncertainty on all measured V2 values. 



5 



J/l/A Elliptic Flow in Pb-Pb Collisions at y/s^ = 2.76 TeV 



The ALICE Collaboration 



^cm 0.3 



- • Pb-Pb ^ = 2.76 TeV, 2.5 <y< 4.0, 1.5 < p T < 10.0 GeV/c 



0.2- 



0.1 







-0.1- 



> 

CD 

CD 



o 
o 
c 

zs 



3.5 



3- 



(a) 



-I — i — i — i — i — i — i — i — i — i — |- 



-I — i — i — i — i- 



1 



(b) 



50 100 150 200 250 300 350 400 

< N ) 
x part ' 

Fig. 3: (Color online) v 2 (a) and (pj) unc0T (b) of inclusive J/y/ with 1.5 < pj < 10 GeV/c as a function of the 
number of participating nucleons in Pb-Pb collisions at ^Ann = 2.76 TeV. 

Figure [2] shows the transverse momentum dependence of the inclusive J/i/a V2 for semi-central (20%- 
40%) Pb-Pb collisions at -y/sNN = 2.76 TeV. The vertical bars show the statistical uncertainties while the 
boxes indicate the point-to-point uncorrected systematic uncertainties, which include the uncertainties 
from the signal extraction, the v^ 8 shape and from the reconstruction efficiency. The global correlated 
relative systematic uncertainty on the event plane resolution is 1.3%. A non-zero V2 is observed in 
the intermediate transverse momentum range 2 < pj < 6 GeV/c. Taking into account statistical and 
systematic uncertainties the combined significance of a non-zero V2 in this pj range is 2.7a. At lower 
and higher transverse momentum the inclusive J/y is compatible with zero within uncertainties. 

To study the centrality dependence of the V2 we select J/i/a with 1.5 < pj < 10 GeV/c for which the 
signal to background ratio as well as the observed V2 are maximized. Since the initial spatial anisotropy 
for head-on collisions is small, the expected elliptic flow is also small. Therefore, we do not consider 
the 0%-5% centrality range. Figure [3] (a) shows V2 for inclusive J/i/^with 1.5 < pj < 10 GeV/c as a 
function of the number of participating nucleons in Pb-Pb collisions at ^Ann = 2.76 TeV. The average 
number of participant nucleons (N pait } for the centrality classes used in this analysis are derived from 
a Glauber model calculation II 1 811 . The vertical bars show the statistical uncertainties while the boxes 
indicate the point-to-point uncorrected systematic uncertainties, which in addition to those discussed 
above also include the uncertainty from event plane resolution determination. The measured V2 depends 
on the pj distribution of the reconstructed J/i/a, which could vary with the centrality of the collision. 
Therefore, (px) uncor of the reconstructed J/i//Ts also shown in Fig. [3](b) as a function of (N paTt }. For 
the most central collisions, 5%-20% and 20%-40% the inclusive J/yv 2 for 1.5 < pj < 10 GeV/c are 
0.101 ±0.044(stat.) ±0.032(syst.) and 0.116 ± 0.045 (stat.)± 0.041 (syst), respectively. The combined 
significance of a non-zero V2 measurement is 2.9a. For most peripheral Pb-Pb collisions, i.e. the two 
classes with low (N part ) values, the V2 is consistent with zero within uncertainties. Although there is a 
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Fig. 4: (Color online) Inclusive J/y V2(pt) for semi-central (20%-60%) Pb-Pb collisions at -^/snn = 2.76 TeV. 
The V2 was measured in the pj ranges: 0-2, 2-4, 4-6 and 6-10 GeV/c and the points are located at the measured 



(.Pt) 1 " 10 ™- Calculations from two transport models [22] and [23] in the same kinematic range are also shown (see 
text for details). 

small variation with centrality, the (pj) anco1 stays in the range 3.0-3.3 GeV/c indicating that the bulk 
of the reconstructed J/y are in the same pj range for all centralities. Thus, the observed centrality 
dependence of the V2 for inclusive J/y with 1.5 < pj < 10 GeV/c does not result from any bias in the 
sampled pj distributions. For J/y with pj < 1.5 GeV/c (not shown) the V2 was found to be compatible 
with zero within one standard deviation for the four centrality classes. The (pj) uncor: ranges from about 
0.75 to 0.9 GeV/c. 

To allow a direct comparison with current model calculations, the inclusive J/yv2(pj) was a l so cai ~ 
culated in a broader centrality range, namely 20%-60%. Figure |4] shows the inclusive J/y V2{pj) for 
non-central (20%-60%) Pb-Pb collisions at -^snn = 2-76 TeV. In this broader centrality range, the mea- 
sured \>2 signal in the pj range 2^1 GeV/c deviates from zero by 2a. The same trend of V2(pj) is 
observed in the 20%-60% and in the 20%^10% centrality classes. This trend is different from the STAR 
measurement fl5ll at lower collision energy, which is compatible with zero for pj > 2 GeV/c albeit in 
somewhat different (10%-50% and 0%-80%) centrality ranges. Also shown in Fig. @]are two trans- 
port model calculations that include a J/i(/ regeneration component from deconfined charm quarks in 
the medium 11221 12311 . In both models about 30% of the measured J/y in the 20%-60% centrality range 
are regenerated. On the one hand, thermalized charm quarks in the medium transfer a significant ellip- 
tic flow to regenerated J/y. On the other hand, initial J/y emitted out-of-plane traverse a longer path 
through the medium than those emitted in-plane resulting in a small apparent V2. The predicted maxi- 
mum V2 at pj ~ 2.5 GeV/c results from an interplay between the regeneration component, dominant at 
lower pj, and the initial J/y component which takes over at higher pj. The first model 02211 is shown for 
the hypothesis of thermalization (full line) and non-thermalization (dashed line) of b quarks. The LHCb 
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Collaboration measured the fraction of J/y from B hadron decays in pp collisions at y/s = 2.76 and 7 



TeV M241I25I1 in the rapidity acceptance used for this measurement. At 7 TeV this fraction increases from 
about 7% at p T ~ to 15% at p T ~ 7 GeV/c, while at 2.76 TeV it is about 7% for p T < 12 GeV/c. In 
Pb-Pb collisions this fraction could increase to a maximum of 11% if the B hadron R^a is equal to 1. 
If b quarks do thermalize then their elliptic flow will be transferred to B mesons at hadronization and 



to the J/y/ at the B meson decay. In the second model 112311 (dash-dotted line) only the case assuming 
thermalization of the b quark is shown. Both models are able to qualitatively describe the pj dependence 
of the V2 and the nuclear modification factor of inclusive J/i/a [4|]. 

In summary, we reported the ALICE measurement of inclusive J/i/a elliptic flow in the range < pj < 10 
GeV/c at forward rapidity in Pb-Pb collisions at ^/s NN =2.76 TeV. For semi-central collisions indica- 
tions of a non-zero J/i/a V2 are observed in the intermediate pj range. This measurement complements 
the results on the J/i/a nuclear modification factor, where a smaller suppression was seen at low trans- 
verse momentum at the LHC compared to RHIC. These results suggest that a significant fraction of the 
observed J/y/ is produced from deconfined charm quarks in the QGP phase. 
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